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ABSTRACT 

We report the result of a search for the infrared counterparts of 37 plane- 
tary nebulae (PNs) and PN candidates in the Spitzer Galactic Legacy Infrared 
Mid-Plane Survey Extraordinaire II (GLIMPSE II) survey. The photometry 
and images of these PNs at 3.6, 4.5, 5.8, 8.0, and 24 /zm, taken through the In- 
frared Array Camera (IRAC) and the Multiband Imaging Photometer for Spitzer 
(MIPS), are presented. Most of these nebulae are very red and compact in 
the IRAC bands, and are found to be bright and extended in the 24 /an band. 
The infrared morphology of these objects are compared with Ha images of the 
Macquarie-AAO-Strasbourg (MASH) and MASH II PNs. The implications for 
morphological difference in different wavelengths are discussed. The IRAC data 
allow us to differentiate between PNs and H II regions and be able to reject 
non-PNs from the optical catalogue (e.g., PNG 352.1— 00.0). Spectral energy dis- 
tributions (SEDs) are constructed by combing the IRAC and MIPS data with 
existing near-, mid-, and far-IR photometry measurements. The anomalous col- 
ors of some objects allow us to infer the presence of aromatic emission bands. 
These multi-wavelength data provide useful insights into the nature of different 
nebular components contributing to the infrared emission of PNs. 

Subject headings: infrared: ISM — planetary nebulae: general — stars: AGB 
and post-AGB 



Figures 2-43 are available on the ApJ website or on request from the authors. 
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INTRODUCTION 



The traditional view of planetary nebulae (PNs) is that it is a shell of gas ionized by a 
hot central star and exhibits an emission-line spectrum. The theoretical understanding that 
PNs descend from asymptotic giant branch (AGB) stars and should contain remnants of the 
dust envelope led to the realization that PNs should also be bright infrared objects (IKwok 
19821 ). This was confirmed after the Infrared Astrono mical Satellite (IRAS) all-sky survey 
found that over 1000 PNs show far infrared emission (IPottasch et al.lll984t IZhang fc Kwok 
199ll ). The typical color temperature of the dust component in PNs is between 100 and 200 
K, warmer than the dust component in H II regions but cooler than the circumstellar dust 
envelopes of AGB stars. This is consistent with the fac t that the d ust envelopes in PNs 
represent the dispersing remnants of the AGB envelopes (lKwoklll990l ). Infrared emission is 
now recognized as a defining observational properties of PNs. 

However, the location of the infrared emitting dust component remains poorly known. 
Only recently has the angular resolution of mid-infrared cameras on large ground-based 
telescopes with adaptive optics become high enough to image the distribution of dust in 
PNs. Space observations therefore offer a better opportunity to survey a large number of 
PNs. The Infrared A rray Camera (IR AQ on the Spitzer Space Telescope has an angular 
resolution < 2 arcsec (IFazio et al.ll2004l ) a nd its capab i lity in imaging the dust component 
of PNs has been clearly demonstrated by iHora et al.l (12004 ) . Since the peak of the dust 
emission from PNs occurs at wavele ngths > 20 /im, the Multiband Imaging Photometer 
for Spitzer (MIPS; iRieke et al.l 12004 ). although offering less angular resolution than IRAC, 
would provide useful information on the dist ribution o f the m ain cold dust component. Based 
on the MIPS observations at 24 and 70 /im, ISu et al.l (120071 ) found possible evidence for the 
presence of a debris disk around the central star of the H elix Nebula. Other exam ples of 
PNs studied by Spitzer includes NGC2346 JSu et al. I2004T). NGC650 dUeta Ibood ). and a 
sample of PNs in the Large Magellanic Cloud (LMC; IHora et al.l 120081 ) . 



One of the advantages of IR observations is that they are hardly affected by dust red- 
dening and are therefor e ideal to study PN s heavily obscured by interstellar dust. From 
the IRAC observations, I Cohen et al.l (j2005al ) identified an extremely red PN G313.3+00.3 
which is optically invisible. This raises the possibility that there are many obscured PNs 
in the Galaxy, for which Spitzer IR observations could explore. A high-spatial-resolution 
and high-sensitivity Ha survey of the southern Galactic plane has been undertaken using 
the Anglo- Australiam Observa tory UK Schmidt Telescope (AAO /UKST) to search for very 
faint PNs (jParker et al.l 119991 ) . This has led to the discovery of over 1000 new PNs and 
PN candidates (which we refer to hereafter simply as PNs). The results of this survey was 
published as the Macquarie/AAO/Strasbourg Ha Planetary Galactic Catalogue (MASH) 
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f lParker et al.l 120061 ) and the MASH II supplement (IMiszalski et al. 1120081 ). These new PNs 
are typically more obscured than known Galactic PNs and warrant further investigation at 
IR wavelengths. 

The Spitzer GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire; 
see Benjamin et al. 2003 and Churchwell et al. 2009) survey provides IRAC images with 
high resolution and sensitivity at wavelengths from 3.6 to 8.0 /iin, allowing us to perform 
high-quality infrared imaging and photometry of PNs. The first-release GLIMPSE I data 
cover a Galactic region of \l\ = 10°-65° and |b| < 1° (total area ~ 220 square degrees). 
Com prehensive GLIMPSE I-based studies of PNs have b een presented by some author s 



Cohen et al. 



2007bl : IPhillips fc Ramos-Lariosll2008al lbl: iRamos-Larios fc Phillip^ l2008h . 



Kwok et all fl2008l . hereafter Paper I) imaged 30 PNa^l detected in the GLIMPSE I survey 
and constructed IR spectral energy distributions (SEDs) which clearly show a distinction 
between the photospheric, nebular bound-free, and dust emission components. 

In this paper, we extend the work of Paper I to the GLIMPSE II survey area that covers 
the Galactic bulge regions. We present the photometry and images of 37 PNs detected in 
the GLIMPSE II survey. These results are compared with those of the GLIMPSE I PNs. 
To better explore nebular dust emission, we also take into account the 24 fim data from 
the Spitzer Legacy program MIPS Inner Galactic Plane Survey (MIPSGAL). This paper is 
organized as follows: in Section 2 we briefly introduce the GLIMPSE II observations and the 
Spitzer data used in this study; in Section 3 we describe the data processing; in Section 4 
we present the Mid-IR images and our measurement results of fluxes, sizes, and SEDs for 
individual PNs; the astrophysical implications of our results are discussed in Section 5; and 
Section 6 summarizes the main conclusions. 



2. OBSERVATION 

The data used in this paper were taken from the GLIMPSE II Version 2.0 Data Release!. 
The GLIMPSE II survey covers longitudes within 10° of the center regions of the Galaxy 
(|6| < 1° for \l\ = 10°-5°, \b\ < 1.5° for \l\ = 5°-2°, and \b\ < 2° for \l\ = 2°-0°; total 
area ~ 54 square degrees) at 3.6, 4.5, 5.8, and 8.0 /xm with resolutions from 1.6"-1.8", but 
excludes an area of ~ 3 square degrees at the Galactic center observed by the PID 3677 GO 
program (which however has been included in the GLIMPSE II v2.0 mosaics). The survey 



2 Note that some objects in Paper I are actually compact H II regions, but are referred to as 'PNs' for 
convenience of denotation (their non-PN nature has been discussed in text). 
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http://www.astro.wisc.edu/glimpse/glimpsedata.html 
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regions are roughly shown in Figured! The observations were made with three 1.2 second 
exposures at each position during the periods between September 2005 and April 2006. The 
procedures of data reduction are similar to those described in Paper I. The GLIMPSE II 
data products include mosaic images, a highly reliable Point Source Catalog (GLMIIC), and 
a more complet e Point Source Ar c hive ( GLM IIA). Further details of the GLIMPSE II data 
can be found in IChurchwell et al.l (120071 ) and iMeade et al. I (120081 ) . 

We also made use of the data from the MIPSGAL surve 

jB The MIPSGAL observa- 
tions were carried out at several epochs between 2005-2006 utilizing the MIPS instrument 
aboard Spitzer, aiming at providing a longer-wavelength complement to the GLIMPSE sur- 
vey. An area of 278 square degrees in the inner Galactic plane, completely overlapping the 
GLIMPSE II observations, was imaged at the 24 and 70 /im bands. In this study, we only 
used the publicly available 24 fim images of the same region as GLIMPSE II. These images 
have a basic sensitivity of ~ 110 /iJy and a field of view (FOV) of ~ 5.4' x 5.4'. The detector 
for MIPS24 is 128 x 128 pixel Si:As arrays with a bandwidth of 4.7/iin. The 24 /im mosaics 
have a saturation limit of 1700MJy/sr for extended source and a resolution of 6". The data 
calibration was performed by both the Spitzer Science Center (SSC) and the MIPSGAL 
team. A goo d summary of the current status of the MIPS24 data sets can be found in 



Cohen! (120091) . The data re ported in this paper were taken from the MIPSGAL v3.0 Data 



Delivery (jCarey et al.l 120081 ) . which provides the 24 /im mosaics (1.1° x 1.1° in size) for all 
the survey regions. 



DATA PROCESSING 



We compared the GLIMPSE II data sets with the Edinburgh / AAO /Strasbourg Cata- 
logue of Galactic Planetary Nebulae (version 1.0; |Parker et alkoOlh and the MASH/MASH II 



catalogue (hereafter referred to as MASH PNs). These optical catalogues include the new 
PNs identified in the AAO/Strasbourg Ha survey of the Milky Way and a small percentage 
of previous known Galactic PNs. We first selected a subsample of 71 MASH PNs lying in 
the GLIMPSE II survey region, among which 10 were taken from the MASH II catalogue. 
We then undertook careful visual examination of the GLIMPSE II images in all the four 
IRAC bands using the DS9 software paekagdj If a FN is visible at at least one of the bands, 
we extracted a portion of the four-band images centered on the coordinates of this PN and 
measured its integrated fluxes. 



4 see http://mipsgal.ipac.caltech.edu/ 

5 DS9 is developed by Smithsonian Astrophysical Observatory. 
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The total fluxes of these PNs were obtained using the same method as in Paper I. We 
first deduced the total on- and off-source fluxes using two apertures which have an identical 
size and were respectively centered on the PNs and in a position near the sources. Then 
the sum of all the fluxes of sources in the point-source archive within each aperture was 
obtained and subtracted from the on- and off-source fluxes respectively. The PN fluxes were 
finally determined from the on-source net fluxes subtracted off-source net fluxes. We have 
taken care to avoid over- and under-subtraction of the fluxes of point sources. For saturated 
sources, we only estimated lower limits of the fluxes. Since most of the GLIMPSE II PNs 
have well-defined boundaries, we also measured the major and minor diameters of each PN. 
For a few sources which have diffuse structures, the sizes of the bright parts were roughly 
estimated. 

The same procedure is also applied to the MIPS24 data. There are several MASH PNs 
which have prominent 24 /im emission but have no IRAC counterparts. We do not include 
these PNs in this analysis. Given the fact that the MIPS24 images have a resolution about 
three times lower than the IRAC ones, in some cases field stars cannot be resolved from the 
PNs in the MIPS24 images, and thus the flux measurements should be treated with some 
caution. The 24 /im mosaics show some artificial structures around bright sources. We have 
attempted to minimize the effects of artifacts on the measurements of fluxes and sizes. Some 
PNs show faint extend halos which are likely to be real. We did not take into account these 
faint halos in the flux estimation since the fluxes of these halos are negligible compared with 
those of the main nebular regions. 



4. RESULTS 

We find that 37 MASH PNs (including five from the MASH II catalogue) have visible 
IRAC counterparts, among which nine are previously known Galactic PNs (Tabled]). The 
detection rate of the GLIMPSE II PNs is 52%, about three times higher than that of the 
GLIMPSE I PNs (18%; Paper I). This is partly ascribed to the fact that the interstellar 
extinction toward the Galactic center is heavier and thus the optical observations only re- 
vealed those relatively bright PNs within the GLIMPSE II field. Most of the GLIMPSE II 
PNs appear to be compact and spherical, probably suggesting that they are generally more 
distant than GLIMPSE I PNs although we cannot completely rule out the possibility that 
some GLIMPSE II PNs are young and intrinsically compact. All the GLIMPSE II PNs are 
strong 24 /im sources, some of which are even saturated. 

Figure. [1] gives the approximate survey coverage of the GLIMPSE I/II projects overlaid 
with the positions of MASH PNs. A total of 12 GLIMPSE II PNs are located within the 



- 6- 



regions of |6| < 1°. This corresponds to a surface density of iVpN = 3 PNs per 10 deg 2 , a factor 
of six lower than that of GLIMPSE II PNs in the regions of \b\ > 1°, where JV PN = 18 PNs 
per 10 deg 2 . This is due to the strong IR background emission w hich hampers the det ections 
of those PNs close to the Galactic plane, as one can clearly see in lMeade et al. I (120071 ) (their 
Figure. 1). This suggests that IRAC counterparts are more easily found outside the galactic 
plane. 

In Figure. [T] we also plot the coverage of the GLIMPSE 3D survey (PI: Benjamin; PID= 
30507), which images 112 square degrees in a series of regions with higher Galactic latitudes in 
all the four IRAC bands aiming at investigating the vertical stellar and interstellar structure 
of the inner Galaxy. We suppose that more PNs can be discovered in the GLIMPSE 3D 
field and it is interesting to compare the properties of these PNs located at higher Galactic 
latitudes with those of the GLIMPSE I /II PNs. Our work on the PNs detected in the 
GLIMPSE 3D survey is in progress and will be the subject of a forthcoming paper. 

The IRAC images of the 37 PNs are displayed in Figure [214381 where we also overplot 
the contours of Ha emission for comparison. The 24 /xm images of these objects are displayed 
in Figure I3"9TT4"31 In order to better show the features at different bands, we have adjusted 
the imaging depths as well as contrasts to create these color images. It is evident from 
these figures that most of these PNs appear to be redder and more extended than the 
surrounding field stars. The mid-IR (MIR) photometry results are presented in Table. [2j 
The first column lists the source names which are designed as "PNG lb" following the IAU 
convention. Column 2 and 3 give the source right ascensions and declinations (J2000), which 
are mostly taken from the MASH catalogues. In a few rare cases the MASH catalogues give 
incorrect coordinates (e.g. PNG 000.9 + 01.8) the corrected values are given in this table. 
Column 4-8 give the spatially integrated fluxes in the 3.6-24 /xm bands. Due to the very 
bright background emission that tends to overwhelm the near-IR emission from the PNs, 
most of the PNs near the Galactic center were detected by IRAC only in the 8.0 /xm band. 
The measured source sizes are presented in Column 9-13. For comparison, column 14 gives 
the Ha sizes taken from the MASH catalogues. For the PNs that are detected in all the four 
IRAC bands, we calculated the 3.6-8.0 /xm SED slopes using apRAc = d\og(\F\)/d\og(\). 
The resultant aiRAC values are given in column 15. 



4.1. Spectral Energy Distributions 

To better constrain the IR spectral properties of the GLIMPSE II PNs, we have searched 
other data archives, including the Two Micron All Sky Survey (2MASS), Deep Near Infrared 
Survey of the Southern Sky (DENIS), Midcourse Space Experiment (MSX), and the IRAS 
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Point Source Catalogue (PSC). Data from these catalogues together cover the spectral re- 
gion from 0.82 to 100 /im. All the available images have been checked visually to ensure 
that the correct counterpart is used. Table [3] tabulates the magnitudes and fluxes of these 
measurements. As these other surveys all have lower sensitivities than IRAC, in general only 
the brightest sources in the GLIMPSE II sample are detected in the other surveys. Since 
IRAC also offers better angular resolution, all the sources are unresolved in these archives, 
with the only exception of PNG 352.1 — 00.0. As PNG 352.1 — 00.0 has a large angular size, 
the data given in these point source catalogues only represent the lower limits of fluxes. We 
therefore do not include this object in Table [31 Combining our IRAC and MIPS photometry 
with these archive data, we have constructed SEDs for the PNs listed in Table [3] and the 
results are displayed in Figure | 



A glance at Figure HH shows that almost all the GLIMPSEII PNs are red objects with 
a rising spectrum between 1 and 20 /im. Their SEDs are in general consistent with the 
presence of a cold dust component of temperature lower than 300 K. This confirms the IRA S 
observations that PNs are generally surrounded with cold dust shells ( iZhang fc Kwok!ll99ll ). 
Some sources have a separate emission component at short wavelengths (< 3/im), which 
can be attributed to contr ibutions from the ph otospheric, nebular bound-free, or hot dust 
emission components (see, IZhang Kwok!ll99ll . for a detailed discussion of SEDs of PNs). 



The SED slope c^rac can provide quantitative description of the SED shape. We have 
determined «irac through least-squares fit to the four IRAC bands, and their values are 



listed in Table EJ The SED slopes for all the GLI MPSE II PNs are a 



larger than those found for stellar photospheres (ILada et al. 



2006 



ac > — 1, substantiall y 



Chavarria et al.ll2008l ). 



Based on their study of young stellar objects (YSOs), IChavarria et al.l (120081 ) introduced 
a classification scheme: Class I for «irac > (star plus infalling envelope), Class II for 
—2 > Chirac < (star plus a disk), and Class III for c^rac < — 2 (post-T Tauri stars or 
main-sequence photospheres). According to this classification, we found that among the 20 
PNs with Chirac measurements 80 perc ent are categorized as Cla ss I (star with envelope), and 
20 percent as Class II (star with disk). Ide la Fuente et al.l (120081 ) determined the c^rac values 
for 19 ultracompact (UC) H II regions, which are comparable with those for PNs. Although 
further statistical work based on a larger PN sample needs to be performed to determine 
whether the classification of YSO can apply to PNs, we show here that Chirac can serve as a 
useful indicator to investigate nebular properties. Nevertheless, it should be borne in mind 
that the slopes of PNs' SEDs can be significantly affected by various emission features (as 
discussed below) and it is difficult to reach solid conclusions without spectroscopic studies. 

Although the SEDs are useful to distinguish contributions from the photospheric, nebu- 
lar, and dust components, the presence of various emission features complicates the IRAC/MIPS 
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based studies of nebular IR properties. These possible emission features include: the con- 
tributions by the 3.3, 6.2 and 7.7 /im aromatic infrared bands (AIB) to the 3.6, 5.8, and 
8.0 /tin IRAC bands, molecular hydrogen vibrational-rotational bands to all the IRAC bands, 
recombination lines (Pfy 3.7/im, Bra 4.05 /im, and Pfa 7.5 /im) to the 3.6, 4.5, and 8.0/im 
bands, and atomic collisionally excited lines to the 4.5-24 /tm bands (see also Paper I). It is 
impossible to accurately estimate the contributions of the se features without spectroscopy . 
The ISO spectrum of PNG 359.3-00.9 (Hb 5) is available JPottasch & Surendiranathll2007h 
and is plotted in Figure. HH Figure. H5] gives an expanded view of this spectrum overlaid 
with the spectral response curves for the IRAC and MIPS 24 /im bands. This figure shows 
that the contributions from these features to these bands are significantly large to change 
the MIR colors. 



4.2. MIR Colors 



In Figure H6] we plot the [3.6] — [4.5] versus [5.8] — [8.0] color-color diagram for the 
GLIMPSE II PNs. For comparis on, this figure also plots the colors of th e GLIMPSE I PNs re- 
ported in the literature (Paper I; Hora et al. 2004 ; Cohen et al. 2007al lbl; |Phillips Sz Ramos-Larios 



2008 



. There is no obvious color difference between the GLIMPSE I and the GLIMPSE II 



PNs in this diagram. These PNs are very red and clustered around a position near (2, 1). For 
the PNs with MIR emission dominated by contributions from the central stars, the IRAC 
colors are close to the (0, 0) point. We have not attempted to correct for interstellar extinc- 
tion towards these PNs (which is supposed to be small at IRAC wavelengths) but show a 
reddening vector of Ay = 10 mag in this figure. We note that the [5.8] — [8.0] color index de- 
creases with increasing reddening because extinction in the 8.0 /tm band is primarily caused 
by a silicate absorption feature peaking at 9.7 /tm. Although one expects Galactic bulge 
PNs to suffer a higher degree of interstellar extinction, but this is not evident in Figure H6j 
The GLIMPSE II PNs are indistinguishable from the GLIMPSE I PNs in the color-color 
diagram. Furthermore, the data scattering direction in Figure H6] is roughly perpendicular 
to the reddening vector. We thus conclude that MIR extinction is negligible and does not 
affect the colors of the PNs in our sample. The IRAC colors of PNs mainly reflect the 
color temperatures and the relative contributions from different components including the 
photosphere, the ionized nebula, the dust envelop, and various emission features. 

We have calculated the simulated colors using blackbody models convolved with the 
spectral response functions, and the results are plotted in Figure H6l Using the simulated 
blackbody colors as a reference, we can see that the [3.6] — [4.5] colors of PNs suggest 
significantly higher temperatures than those by the [5.8] — [8.0] colors. This is indicative of 
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the substantial contribution of emission from AIB emissions, ionized gas, and/or photosphere 
to the 3.6 nm band. The relative contribution of the latter two have been illustrated in 
Zhang fc Kwokl Jl99lh . 



In spite of the fact that PNs have different colors with stellar objects, the color-color 
diagram alone does not enable to completely distinguish PNs from other MIR sources. From 
Figure H6] we can find that the colors for PNs partly overlap with those for Class I YSOs (i.e. 
star with envelope) and the PNs with strong emission from central stars might be confused 
with Class II YSOs (i.e. star with disk). Moreover, some H II regions are located within the 
same position as PNs in the color diagram. Therefore, despite the color-color diagram can 
serve as a tool to sample PN candidates, detailed investigations of morphologies and spectra 
are essentially required for accurate PN designation. 

Figure 1471 shows the [3.6] — [8.0] versus [8.0] — [24] color-color diagram. The [3.6] — [8.0] 
colors are in the range from 1-5, and the [8.0] — [24] colors are in the range from 3-7. 
The IRAC-MIPS colors for these Galactic PNs are similar to those found for LMC PNs by 
Hora et al.l (120081 . see their Figure 5). For the same reason as in Figure |4"6| the PN colors 
cannot be explained by single-temperature blackbodies. Comparing Figure 1471 with the co lor- 



color diagram for the members of the IC 348 cluster presented by lLada et al. I (120061 . see 
their Figure. 7), we find that the distribution of PNs in this diagram are clearly centered 
at a position redder than IC 348 cluster sources, even for the YSOs with thick disk. This 
suggests that the 24 /im excess is a prominent characteristic of PNs. 

In Fig. HH]we plot the IRAC color-magnitude diagrams ([3.6] versus [3.6] — [8.0] and [8.0] 
versus [3.6] — [8.0]). These diagrams show that the GLIMPSE II PNs generally have lower 
integrated fluxes than the GLIMPSE I PNs, although their colors are similar to each other. 
These diagrams can be compared with those for LMC PNs (IHora et al.ll2008l . Fiugre 6). We 
find that the distribution of Galactic PNs in the color-magnitude diagrams perfectly follows 
those of LMC PNs in shape although the galactic PNs are nearer and hence brighter. 



4.3. Individual Objects 



PNG 00.0—01.3. — This so urce was detected by IRAC in the 5.8 and 8.0 /xm bands only 
(Figure [2]). iParker et ajj (120061 ) classifie d it as a 'Likely' PN . The object was also detected 
by the NRAO VLA Sky Survey (NVSS; ICondon et al.lll998l ). which measured a flux density 
of 4.5 ± 0.8 mJy at 1.4 GHz. In spite of its faintness, its Ha and 8.0 /im images are visually 
similar and both reveal an oval structure. 



PNG 000.1—01.7. — This source was detected by IRAC in the 8.0 /im band only (Fig- 
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ureBl) . Its PN status has been confirmed by spectroscopy (designated as JaSt85; |Jacoby fe Van de Steene 
2004 ). IVan de Steene fe Jacobyl (120011 ) observed this Galactic bulge PN in the radio con- 
tinuum at 3 and 6 cm and estimated a dis tance of ~ 7.2 kp c. The NVSS survey gives a 



flux d ensity of 2.9 ± 0.6 mjy at 1.4 GHz (ICondon et al.lll998f ). Given large [N II] /Ha flux 
ratio (IParker et al.ll2006l ). it might be a low excitation PN. It is extremely faint in the IRAC 
bands but fairly bright in the MIPS24 band. Albeit with an ambiguous boundary, the size 
of 8.0 /im image is roughly consistent with that of Ha. 

PNG 000.3—01.6. — This source was detected by IRAC in the 8.0 /xm ba nd only (Fig- 

ureHD. Its PN status has been confir med by spectroscopy (designated as JaSt86; |Jacoby fe Van de Steene 




Based the radio observation, IVan de Steene fe Jacobyl (120011 ) estimate d a distance of 



6.6 kp c. The NVSS survey gives a flux density of 4.1 ± 0.5 mjy at 1.4 GHz (ICondon et al. 
19981 ). It is very compact in all the MIR and Ha bands. 



PNG 000.5+01.9.— This source wa s detected by IRA C in the 8.0 /im band only (Fig- 
ure ED- The object was first discovered by lKohoutek I (119941 ) (designated a s K 6-7), and its PN 
statu s was subsequently confirmed by s pectroscopy (designated as JaStl7; |Jacoby fe Van de Steene 
2004 ). IVan de Steene fe Jacobyl (120011 ) estimated the distance to be 7.1 kp c. The NVSS sur- 
vey gives a flux density of 6.8 ± 0.9 mjy at 1.4 GHz (ICondon et al.lll998l ). The 8.0/im and 
Ha images are visually similar and show a compact and oval structure. 

PNG 000.9-01.0.— This source was detected by IRAC in the 8.0 //m band only (Fig- 
ure ED- It is classified as a 'True' PN in the MASH II catalogue. The MIR ( 8.0 and 24 /im) 
and Ha bands show a compact and round structure without a sharp boundary. 

PN G 000.9+01.8 . — Th is source was detected by IRAC in the 8.0 /im band only (Fig- 



ure [7]). IParker et al.l (120061 ) classified it as a 'True' PN. The optical image shows an oval 
structure. From the IRAC image, we note that the object seems to be a core surrounded by 
an ambiguous diffuse halo. 

PNG 001.0-01.9.— This source was de tected by IRAC i n the 8.0 //m band only (Fig- 
ure ED- The object was firs t discovered b y iKohoutek I (120021 ) (designated as K 6-35) and 
classified as a 'True' PN bvlParker et al.l (120061) . The NVSS survey gives a flux density of 
4.8±0.6mJy at 1.4 GHz ( Condon et al. 19981 ). The Ha image shows that it is a bright, com- 
pact, and slightly oval PN with faint outer halo. Its 8 /im counterpart is nearly overwhelmed 
by bright background emission. 

PNG 001.0+01.9. — This source was only marginally detected by I RAC in the 8.0 /im 
band and was completely invisible in the other IRAC bands (Figure [9]). lExter et al.l (120041 ) 



assified it as a 



determined chemical abundances of this object with optical spectroscopy and c 
Type I PN. The NVSS survey gives a flux density of 20.6 + 1.3 mjy at 1.4 GHz (ICondon et al. 
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19981 ). IZhana (119951 ) deduce d the statistica l dista nce to be 3.46 kpc. In their catalogue of 
narrow band images of PNs, ISchwarz et al.l (119921 ) presented Ha and [0 III] images of this 
PN, which exhibit an oval appearance with a size of 74" and of 43", respectively. The 8.0 /im 
image is too faint to obtain a reliable flux. This PN is different from others in this sample 
of PNs in that its emission is significantly more compact at 24 /im than at Ha. 

PNG 001.6—01.1. — This source was detected by IRAC in the 8.0 /im ban d only (Fig- 



urefTUl). Its PN status has been confirm ed by spectroscopy (designated as JaSt97; |Jacoby fc Van de Steene 
20041 ) . IVan de Steene fc Jacobyl (120011 ) estimated the distance to be 5.9 kpc. The 8.0 /im and 
Ha images have a similar size and both exhibit a compact oval appearance. 



PNG 002.0+00.7.— This source was detected by IRAC in all th e four bands (Figure ED. 
The NVSS survey gives a flux density of 4.1 ± 0.5 mJy at 1.4 GHz JCondon et al.lll998l ). It 
is classified as a 'True' PN in the MASH II catalogue. The PN appears to be a point source 
in all the bands. Its IRAC color is only slightly redder than stars. The SED indicates the 
presence of cold dust with T < 300 K (Figure 144)) . 

PN G 002.1-00.9 . — Th is source was detected by IRAC in all the four bands (Fig- 



ur e [T2l) . iParker et al.l (120061 ) classified it as a 'Likely' PN. It was also independently found 
by iKohoutek I (120021 ) . This PN is very bright in all the bands, and exhibits a compact and 
round appearance. The SED indicates a cold dust component with T < 300 K (Figure | 



PNG 002.1—01.1. — This source was detected by IRAC in the 5.8 and 8.0 /im bands 
only (Figure [13]). It is classified as a 'Ture' PN i n the MASH II catal ogue. The NVSS survey 
gives a flux density of 2.1 ± 0.5 mJy at 1.4 GHz (j Condon et al.lll998l ). Its Ha image shows a 
seemingly bipolar structure, which however is hard to tell from its fuzzy 8 /im image. 

PNG 002.2+00.5.— This source was detected by IRAC in all t he four bands (Figu re ED . 
The NVSS survey gives a flux density of 15.2 + 0. 7mJy at 1.4 GHz JCondon et alJll998h . The 
Ha image exhibits an elliptical appearance with resolved internal structure. The IRAC data 
show that it is a very red PN and has a morphology visibly similar to that of the archetypal 
PN NGC 7027, i.e., a limb-brightened, bi-conical shell with enhanced emission at the equator 
(see 



Cox et aU 12002 . for the details of the image of NGC 7027). 



PNG 002.2-01.2.— This source was detecte d by IRAC in the 4 5, 5.8, and 8.0 //m bands 
only (Figure 115]) . It is classifie d as a 'Ture' PN by Parker et al.l ( 120061 ). Strong [N II] and weak 
[O III] lines were detected by lParker et al.l (120061 ). suggesting that it is a low excitatio n PN. 
The NVSS survey gives a flux density of 2.1 ± 0.5 mJy at 1.4 GHz (jCondon et al.lll998l ). The 
8.0 /im and Ha appearances are similar to each other and show a compact and rectangular 
shaped structure. 
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PNG 003.4+01.4.— This source is almost invisible in the IRAC bands (Figure US]). We 
only find an extremely faint counterpart in the 8.0 /im band. Its 24 /iin counterpart appears 
to be a poi nt source and is fai nter compared to other GLIMPSE II PNs. The object is 
described by lParker et al.l (120061 ) as having faint, semi-circular arc around faint central star, 
and is classified as a 'likely' PN. Albeit with faintness in optical and MIR bands, it has a 
sufficient ly strong radio flux . The NVSS survey gives a flux density of 15.2 ± 0.7mJy at 
1.4 GHz (ICondon et al.lll998l ). Its faint 8.0 /im image does not allow us to say much about 
its infrared structure. 

PNG 00 3 .5—01 .2. — This source was detected by IRAC in all the four bands (Figure IT7"|) . 



Parker et al.l (120061 ) detected strong [N II] and weak [O III] lines in the object and classified 
it as a 'Possible' PN, possibly a symbiotic star containing a Mira variable. The object 
has an IRAS counterpart, IRAS 17554 — 262 8. The NVSS survey gives a flux density of 
13.7 + 0.6 mJy at 1.4 GHz (ICondon et al.lll998l ). The IRAC and 24 /im images exhibit a very 
red and bright point source. The SED implies the presence of dust with a temperature of 
~ 300 K (Figure HI]). Moreover, the increase of flux at ~ 100 /im points to the existence of 
another dust component with even lower temperature. 



PNG 003.5+01.3. — This source was detected by IRAC in all the four bands (Figure [T8|) . 
It is classified as a 'True' PN in the MASH II catalog ue. The NVSS survey gives a flux density 
of 24.9 ± 0.9 mJy at 1.4 GHz (ICondon et al.lll998l ). Its IRAC images appear to be a point 
source. The SED shows a marked increase of flux from 6-12 /im, suggesting the presence of 
cold dust of T < 300 K (Figure | 



PNG 003.6— 01.3. — This s o urce w as detected by IRAC in all the four bands (Figure [T9]) . 



It is described by iParker et al.l (120061 ) as a small, bright, round, and 'True' PN with strong 
[O III] and Ha emission. The IRAC images exhibit a similar appearance with Ha image 
although the former seem to be more compact. The SED shows an emission peak around 
1.5 /im (Figure [HJ, probably pointing to the substantial contribution from photospheric and 
nebular bound-free emission. 

PNG 003.7+00.5. — This source is badly blended with a very bright star (Figure 12171) 
which saturates IRAC. Nevertheless, we can still find evidence sho wing extended emiss ion 
in the northwest of this bright field star. The object is described by IParker et al.l (120061 ) as 
a very small, isolated slightly oval nebula with asymmetric enhance opposing edges, and is 
classified as a 'Likely' PN. The 24 /im emission is seriously polluted by the bright field star, 
and thus we only give the upper limit of the flux. 



PNG 004.3-01.4.— This s ource was detected b y IRAC in all the four bands (FigureEID- 
It is classified as a 'True' PN by lParker et al.l (120061 ). The NVSS survey gives a flux density of 
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12.7±0.7mJy at 1.4 GHz (ICondon et alJll998h . This object has a IRA S counterpart, IRAS 
17582 — 2553. The MIPS 24 /xm flux of 1.7 Jy is in excellen t agree ment with the 25 /xm flux 
(1.9 Jy) given in the IRAS PSC. Kistiakowsky fc Helfand I (119951 ) imaged this PN through 
narrowband filter centered at [S III] A9532, which reveals a size of ~ 2". The IRAC and Ha 
images display a compact and roughly round structure. The SED shows a steep increase of 
flux at 10 /xm (Figure 1441). and suggests that the PN might have a dust component with a 
color temperature of ~ 300 K. 

PNG 004.8-01.1.— This source was detected by IRAC i n all the four bands (Fig- 
ure [221) . It is classified as a 'True' PN by iParker et al.l (120061 ). This object has a IRAS 
counterpart, IRAS 17581 — 252 2. The [S III] A9532 narrowband image suggests a size of 4" 
(IKistiakowsky &: Helfand Ill995l ). The smaller size at [S III] compared to those at MIR and 
Ha bands is due to the fact that the former only traces the inner hi gh-ionization regions . 
The NVSS survey gives a flux density of 16.2 ± 3.1mJy at 1.4 GHz (ICondon et al.lll998l ). 
The 8 jum and Ha images show an oval appearance. The SED suggests the presence of cold 
dust with T < 150 K (Figure 144]). 



PNG 006.1+00.8. — This so urce was detec t ed by IRAC in all the four bands (Figure[ 
It is classified as a 'True' PN by IParker et al.l (120061 ) . The 8 /xm and Ha images are similar 
and show a compact and round appearance. 



PNG 352.1—00.0. — This source was detected by IRAC in all the four bands (Figure 
The Ha image shows an arcuate nebula around a star. Much more extended structures are 
revealed by the IRAC images which show a sharp edge to the east and a diffuse tail trailing 
toward the west. Based on its morphology, it is likely to be an H II region and not a 
PN. The cometary structure shown by the MIR images is commonly seen in compact and 
ultracompact H II regions and implies the interaction between the nebula and the ISM. 
We also find that the object is surrounded by very extended (~ lOarcmin across) fainter 
filaments. The 24 /xm emission is dominated by the central bright region and appears as 
a seemingly oval structure. This object is an example showing that these MIR data can 
help for identifying/rejecting PNs. We note that there are two bright IR so urces, IRAS 



1722 1 - 3533 and 2MASS J17253388 - 3536005 (the latter might be a YSO; iFelli et al. 
20021 ) . lying within the nebulosity. 



PNG 352.8—00.5. — This source was detected by IRAC in all the four bands (Fig- 
ure [25]). It has a IRAS counterpart, IRAS 17262 — 3511 and is classified as a 'True' PN 
in the M ASH II catalogue. The NVSS survey gives a flux density of 72.8 ± 2.3 mJy at 
1.4 GHz ([Condon et al.lll998l ). We note that the PN lies within an overlap region between 
the GLIMPSE I and GLIMPSE II surveys and its IRAC images have already investigated by 
Ramos-Larios Phillips! (120081 ) . We refer the reader to their paper for further details of this 
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PN. Its 24 /mi emission is very strong and saturates the detector. The SED shows a smooth 
increase of flux from 1-50 /im, and suggests the presence of cold dust with T = 60-100 K 
(Figure! 



PNG 353.9+00.0. — This source was detected by IRAC in all the four bands (Figure 126]) . 
Optical spectroscopy sugge sts a low excitation class and the object is classified as a 'Possible' 
PN by iParker et al.l (120061 ). The Ha and MIR images appear as a point source. Unlike the 
other PNs, this source has a IRAC color only slightly redder than stars, and has a low 
SED slope ( Chirac = —0.91). The SED indicates to a color temperature of about 1000 K 
(Figure |4*4"|) . or even higher if the interstellar extinction is non-negligible. This value is 
significantly higher than in normal PNs. Consequently, this object is more likely to be an 
emission-line star rather than a PN. 

PNG 355.6—01.4. — This source is almost overwhelmed by the bright background emis- 
sion and was detected by I RAC in the 8.0 /im band only (Figure [27]). It is classified as a 
'True' P N bvlParker et al.l (120061) . The NVSS survey gives a flux density of 5.5 ± 0.6 mJy at 
1.4 GHz (ICondon et al.lll998l ). The 8 /im image is visually consistent with the Ha appearance 
that shows a small bipolar-like structure. 

PNG 355.6+01.4. — This sourc e was detecte d by IR AC in all the four bands (Figure 128]) . 
It is classified as a 'Likely' PN by IParker et al.l (120061 ). The spectroscopy suggests that it 
probably is a very low excitation PN. The NVSS survey gives a flux density of 8.2 ± 0.6 mJy 
at 1.4 GHz (ICondon et al.lll998l ). The Ha and MIR images appear as a point source. The 
SED implies the presence of cold dust with T < 100 K (Figure | 

PNG 356.0 
only (Figure 



-01.4. — This source was detected by IRAC in the 5 .6 and 8.0 /zm bands 
It is classified as a 'True' PN b y IParker et al. (120061) . The NVSS survey 

The Ha and IRAC 



gives a flux density of 4.8 ± 0.5 mJy at 1.4 GHz (j Condon et al 
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images show that it is a faint, compact, and elliptical PN. 

PNG 356.9+00.9.— This so urce was detected by IRAC in all the four bands (Figured)]). 
It is classified 'True' PN bvlParker et al.l J2006J). The NVSS survey gives a flux density 
of 12.5 + 0.6 mJy at 1.4 GHz (jCondon et al.lll998l ). The Ha image appears as a point source. 
However, extended emission is clearly resolved in the IRAC images which exhibit a round 
structure. The PN is very red and bright in the MIR bands. The SED implies the presence 
of cold dust with T < 150 K (Figure [H 



PNG 357.4—01.3. — This source was detected by IRAC in all the fo ur bands (Figu re I3~TP . 
Given its large size (~ 5arcmin), the object is unique in our sample. IParker et al.l (1200 ll ) 
suggested that it might be a highly evolved PN. As the extended emission region is very 
diffuse, we cannot obtain reliable integrated flux and can only roughly estimate the size. 
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The optical image s hows that it is a ring nebula around the Wolf-Rayet (WR) star WR 101. 
Cappa et al.l (120021 ) reported the radio image of this nebula at 1465 MHz taken through the 
VLA and obtained an ionized mass of about 230 M . Therefore, this object is unlikely to be 
a PN, and is probably a WR nebula consisting of mass lost as well as interstellar material 
swept-up by the stellar wind from an evolved massive star. We find that the Ha and 24 /im 
images have a high resemblance. However, compared to that in Ha and 24 /im, the ring 
structure is less well defined and fuzzier in the IRAC bands, which also reveal more diffuse 
emission extending to the northwest. The position of the bright clump in the southeast as 
shown by IRAC slightly differs from that of the Ha counterpart which appears to be closer 
to the ionized source. This suggests that the nebulosity is optically thick to the ionizing 
photons, and probably the dust grains close to the ionized source have been destructed by 
photodissociation of the central star. 

PNG 357.5+01.3. — This sourc e was detecte d by IR AC in all the four bands (Figure 132]) . 
It is classified as a 'Likely' PN by iParker et al.l (120061 ). The Ha image appears as a point 
source. Extended emission is clearly resolved in the IRAC images, which show that it is a 
red and round nebula. 

PNG 357.7+01.4. — This source was detected 



jy IR AC in all the four bands (Figure [33 



It is classified as a 'True' PN by IParker et al.l (120061 ). The Ha image shows a slightly 
elongated structure. The IRAC images only reveal a bright core and the object is the most 
compact PN in our sample. 

PNG 358.2—01.1. — This source was detected by IRAC in all the four bands (Figure[3 



This PN has been discovered before (e.g. ICahn fc Kalenll97ll; lAcker et al.lll991l) . The NVSS 



be 144,000 K. 



Maciel 



survey gives a fl ux density of < 2.5 mJy at 1. 4 GHz (jCondon et al.lll998l ). Using the energy- 
balance meth od. iPreite-Martinez et al.l (119891 ) obtained the temperature of the central star to 



19841 ) determined a distance toward this PN of 1.1 kpc. It is described 



by IParker et al.l ( 20011 ) as a bright, compact, possibly bipolar PN with arc features. The 
IRAC images have a visually similar appearance with the Ha image and show arc features 
and fuzzy extended structure. 



PNG 358.8—00.0. — This source was detected by IRAC in all t he four bands (Figur e! 
The NVSS survey gives a flux density of < 2.5 mJy at 1.4 GHz (jCondon et al.lll998l ). Al- 
though the object has a IRAS counterpart , IRA S 17395 — 2950, we cannot find the flux 
measurements in the IRAS PSC. IVolk et al.l (Il99ll ) presented the low-resolution IRAS spec- 
trum and classified the PN into Group I which represents noisy or incomplete spectra. The 
Ha image shows a small and fuzzy nebula. It has a large ai^Ac value, and is redder than all 
the other GLIMPSE II PNs. The MIR images are remarkably extended and reveal intriguing 
structure. The bright region has a C-shape. We find that it has a large dust envelope with 
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a cavity centered on the central star and diffuse emission extending to the west. Thus it 
is unlikely to be a PN. The SED implies the presence of cold dust with T < 150 K (Fig- 
ure SI]). From the IRAC images, we can also see a bipolar nebulosity in the northwest to 
PNG 358.8—00.0. This object has no visibl e optical counterpar t and is probably associated 
with a compact radio source discovered by iGray et al. I (119931 ). From the morphology, we 
conjecture that it is likely to be a PN or compact H II region highly obscured by dust ex- 
tinction. Follow-up observations of this GLIMPSE-survey-discovered object are required to 
unveil its nature. 



PNG 359.1-01.7.— This source was detected by IRAC in all t 
This is a previously known PN (e.g. iHenize et al.lll967l ; IWebster 



re fou r bands (Figure 136]) . 
A detailed optical 



1975) 



spectrophotometric study has been recently presented by I Wang &: Liu I (120071 ). The deduced 
chemical abundances suggest that it a Type I P N. The NVSS s urvey gives a flux density 
of 91.9 ± 2.8 mJy at 1.4 GHz JCondon et a l.l ll998h. Izhangj Jl995T) derive d a distance to this 
PN of 3.30 kpc. Using a dynamical method, iGesicki fc Zijlstra I (120071 ) obtained the mass 
of the central star of ~ 0.627 M^. The IRAC and Ha images show a compact, bright, and 



round appearance. However, IGesicki et al. I (120031 ) found that this object has complicate 



line profiles of H I, [N II], and [O III], and suggested the presence a bipolar outflow. In the 
IRAC images, we do note that the brightest part is not located in the center of the nebula. 
The SED indicates to a color temperature T < 300 K (Figure | 



PNG 359.2+01.2.— This source was detected by IRAC in all the four bands (Fig- 
ure I3TD. It has been know n as a typical bipolar PN (e.g 



Kohoutek 


1982; 


Ratag et al. 


1990 



Corradi fc Schwarzl Il995l ) and has been suggested as a symbiotic sy stem (ICorradi 



The NVSS survey gives a flux density of 13.5 ± 0.6 mJy at 1.4 GHz (ICondon et al. 



1995) 
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It has a IRAS counterpart, IRAS 17358 — 2854. The Ha image shows a bright core with a 



highly collimated bipolar outflow. The compact core has also been detected by iLee et al. 



( 120071 ) at radio wavelengths. The IRAC images show a less extreme bipolar shape with a 
lower major/minor axis ratio compared to the Ha image. We speculate that the PN has a 
central dust torus therefore enhancing the IRAC emission near the center. Interestingly, the 
appearance revealed by the 24 /mi image is a bipolar nebulosity surrounded by an extended 
halo with a size of 30" x 42". The SED indicates to the presence of cold dust with T < 100 K 
(Figure! 



PNG 359.3—00.9. — This source was detected by IRAC in all the four bands (Fig- 
ure [38]). _TJiis_js_^_well4<nown^ been extensively studied (see, 
e.g. 



Pottasch fc Surendiranath 2007; Montez et al. 2009. for recent literature summaries of 



this object). Optical images show that this PN is composed of two closed-end bipolar lobes, 
a bright compact core, concentric rings and some filaments. It is a high-excitation Type I 
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PN and suffers from heavy extinctio n. The distance to t his PN is not accurately known, 
with estimates ranging from 1-7 kpc (IMontez et al. II2009I ). It is one of the four PNs that 



have hitherto been detecte d with X-ray emissi on, indicating the presence of shock-heated 



gas in the bipolar nebula (IMontez et al. 



179.5 ± 5.4 mJy at 1.4 GHz (ICondon et al. 



20091) . The NVSS survey gives a flux density of 



19981 ) . The IRAC images exhibit the bright com 



pact core, which is so luminous that the 8 /zm band is saturated. Because of the brightness of 
the central core, the bipolar lobes are not obvious from the IRAC images. It is clear that in 
the IRAC wavelengths, the core-to-lobe flux-density ratio is larger than that at Ha. We sug- 
gest that the 8 jum emission is dominated by the dust thermal emission from a central torus. 
The 24 jum band, however, shows more extended emission even when deconvolved using a 
point source function. The 24-/im image suggests an extended cold dust envelope, corre- 
sponding to the ~ 100 K dust component seen in the SED (Figure 1441) . The ISO spectrum 
of Hb 5 is plotted on the SED, showing that ISO flux measurements are in good agreement 
with other measurements. We note that there are two strong high-ionization lines [Ne V] 
24 fim and [O IV] 26 /zm lying in the 24 /xm band, which contribute about 10 percent of the 
total 24 fim flux and are partly responsible for the saturation of the central region in the 
24 /xm image. 



4.4. Nondetections 

The following MASH PNs have no obvious IRAC counterparts: PNG 000.0-01.8, PNG 
000.1 + 01.9, PNG 000.6-01.4, PNG 001.0-01.4, PNG 001.1-01.2 PNG 001.2 + 01.3, PNG 
001.5-01.6, PNG 001.6 + 00.1, PNG 001.6-00.6, PNG 001.6 + 01.6, PNG 001.9 + 01.9, PNG 
002.1 + 01.2, PNG 002.4 + 01.1, PNG 002.5 + 01.3, PNG 002.7-01.4, PNG 004.0-00.4, PNG 
004.8-00.5, PNG 005.5-00.8, PNG 007.2 + 00.0, PNG 009.9-00.2, PNG 354.0-00.8, PNG 
354.3 + 00.5, PNG 354.8-00.5, PNG 355.5-01.1, PNG 355.9 + 00.7, PNG 357.3 + 01.3, PNG 
357.6 + 01.0, PNG 359.1-00.7, PNG 359.2-01.2, PNG 359.2 + 01.3, PNG 359.3 + 01.4, PNG 
359.5-01.8, PNG 359.9-01.8, PNG 359.9 + 01.8. Inmost of the cases, the detection of these 
PNs is hampered by the bright background near the Galactic center. For several objects (e.g. 
PNG 001.1 - 01.2 and PNG 002.7 - 01.4), the IRAC bands exhibit apparent point sources 
near the corresponding positions, but without extended emission around these sources. They 
are likely to be field stars although some of them have a rising flux distribution from 3.8 to 
8.0 /im. An interesting object is PNG 354.8 — 00.5, which is a radio and maser source. Its 
Ha image reveals a faint, circular nebula with a size of 42.0" x 37.0". However, the IRAC 
images show an extremely luminous star, which is bright enough to saturate the detector 
and makes the nebula 'invisible'. The object probably repre sent a evolutionary s tage of PN 



formation around a very evolved AGB star, as suggested by I Cohen et al.l (l2005b[ ). 
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DISCUSSION 



All the objects studied in our sample can be categorized as ionized nebulae with cold 
dust emission. The ionized nature of the objects is confirmed by Ha and radio continuum 
emission. The peak of the dust continuum mostly lies between 10 and 30 /im, suggesting dust 
temperatures of 100-300 K. Since the infrar ed spectra of PNs are rich in ionic emission lines as 
well as solid-state emission bands (see, e.g.. lStanghellini et al.ll2007l ). accurate determination 
of dust temperatures of PNs is very difficult without spectra and our estimates represent 
only rough approximations. 



5.1. Implications for Source Sizes at Different Wavelengths 

It is instructive to compare the morphologies observed in different bands since they could 
trace the spatial distributions of different nebular components. Figure H9] shows histograms 
of source sizes in MIR wavelengths relative to Ha, for which we have deconvolved the point 
spread functions (PSF) simply through 9 S = a/ #obs ~ ^iL (note that this is only a rough 
estimate as the PSFs are non- Gaussian), where 8 S is the intrinsic source sizes, o b s is the 
observed sizes, and 9 Tes is the angular resolution. An inspection of this figure shows that (1) 
the Ha sizes are more similar to those observed in 5.8 and 8.0 /im compared to the other 
MIR sizes, (2) the 24 /im sizes are generally larger than those in the other wavelengths, and 
(3) The 3.6 and 4.5 /im sizes are quite similar to each other and are more compact than those 
in the other bands. 

The optical and 8.0 /im morphologies of a certain PN can be remarkably different. For 
some PNs, the 8.0 /im emission has a larger angular extent than the corresponding optical 
image. This is indicative of the presence of extended dust envelopes which obscure the 
optical emission and emit thermal continuum and AIB features contributing to the 8.0 /im 
band. This can also in part be attributed to the H 2 rotational- vibrational transitions excited 
by shocks. On the other hand, some PNs are more compact at 8 /im than at Ha. An 
important issue is whether this really reflects the different spatial distributions of the different 
components rather than different instrumental sensitivities. In order to investigate this 
problem, in Figure. [50] we compare the Ha and 8.0 /im flux profiles along the major axis of 
the bipolar PN PNG 359.2 + 01.2. While the Ha and 8.0 /im images have comparable angular 
resolutions, it is clear that the 8.0 /zm emission is more centrally concentrated than the Ha 
emission. The most likely interpretation is that the Ha and 8.0 /xm emissions originate from 
different regions, the former from the bipolar lobes and the latter from a nuclear dust torus. 
Examples of infrared- observed equator ial dust torus can be found in the bipolar proto-PN 



IRAS 17441 - 2411 flVolk et al 



20071 ). If we generalize from this example, objects with 
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smaller angular extents in 8.0 fim than in Ha may be bipolar nebulae. 

The 24 fim images of our sample objects are in general well resolved, extended, and 
show a circular appearance with a sharp boundary. It is unlikely that [Ne V] 24.3 fim and 
[O IV] 25.9 /im lines are dominant contributor to the 24 /im flux otherwise the 24 /im image 
would show a more compact appearance than the Ha image. It is more probable that the 
24 /im images correspond to the distribution of cold dust which is the remnant of AGB mass 
l oss. Extended AG B halos have been frequently detected in deep optical images of PNs 



(ICorradi et all 120031 ). If the cold dust grains are plentifully present within these halos, they 



can exhibit strong 24 /im emission. 

For high-excitation PNs, the [Ne V] and [O IV] line emission may significantly enhance 
the 24 fim emission in nebular centers. A recent study of MIPS 24 /im images of Galactic 
PNs (which are generally more extended than objects in our sample) has been reported by 



Chu et all (120091 ) , who find that the sizes and surface brightnesses of Ha and 24 /im emission 



suggest an evolutionary sequence. 

From Figures I3"9"H4"3"| we can see some of our sample objects show faint and extended halos 
and ring structures around their bright cores. While some of these are clearly instrumental 
artifacts, some are likely to be real (e.g., PNG 001.0-01.9, PNG 001.0 + 01.9, PNG 359.1 - 
01.7 etc.) and deserve further study. A recent MIPS24-based study to search for debris 
disk around stars has found 24 /im excess emission in the Helix Nebula and suggested that 
it originates from a dust disk at a distance of 40-100 AU fr om the central star, and that the 



dust is produced by collisions of Kuiper-Belt-like objects (jSu et all 120071 ). Assuming that 
the distances to the GLIMPSE II PNs are typically larger than 1 kpc, we find that the 24 fim 
emission extends to > 5000 AU from the center, more than a factor of 100 times farther 
than the Sun's Kuiper Belt, and thus is unlikely to arise from a Kuiper Belt around low- 
and intermediate-mass AGB stars. However, we cannot completely rule out the possibility 
that Kuiper Belt-like objects may partly contribute to the 24 fim emission in the centers. An 
interesting PN is PNG 001.0 + 01.9 which is the only object showing more compact 24 fim 
emission relative to Ha in our sample. Its 24 fim emission is probably in part associated 
with the central star. 

For the two short-wavelengths IRAC bands (3.6 and 4.5 fim), the contributions from 
the hydrogen recombination line and bound- free emission should not be dominant, because 
if so, the 3.6 and 4.5 fim emission would morphologically match the Ha emission well. The 
spatially compact 3.6 and 4.5 fim emission probably suggest significant contribution from 
photospheric emission. This conclusion is further supported by the SEDs (Figure |4"4"|) and 
color-color diagrams (Figures l4"6l and I4T1) . 
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5.2. IRAC 8.0 /im versus MSX 8.3 /im Integrated Fluxes 



Reach et al.l (120051 ) investigated the absolute calibration of IRAC. They found that the 
extended emission surface brightness is incorrectly calibrated and should be multiplied by 
an effective a perture correction fa ctor (0.737 for the 8.0 /im band). This was subsequently 
established by lCohen et al.l (j2007al ) that the IRAC 8.0 /im surface brightness should be scaled 
down by a factor of 0.74 ± 0.07 based on comparison between the IRAC 8.0 /im and MSX 
8.3 /im fluxes of a sample of 43 H II regions. This study predicts a IRAC8.0/M5X8.3 flux 
ratio of 1.14 ± 0.02 for diffuse H II regions based on the different con tributions from AIB 
emission in the two different bandpasses. A similar study of 14 PNs ( jCohen et al.l l2007bl ) 
also confirm their previous results. Here we present an independent examination of the 
IRAC 8.0 /im calibration using the integrated fluxes of the GLIMPSE II PNs which are very 
compact diffuse sources. 

A total of 10 GLIMPSE II PNs have MSX 8.3 /im detections. Figure M compared 
the spatially integrated fluxes detected in the IRAC 8.0 /im and MSX 8.3 /im bands. For 
comparison, we also overplot the values for GLIMPSE I PNs. A good correlation between 
the two fluxes is clearly shown. The PN with the largest IRAC/ MSX flux discrepancy is 
PNG 003.5 + 01.3, which has a large MSX 8.3 /tm flux compared with the IRAC 8.0 /im 
flux. Combining the GLIMPSE I and GLIMPSE II measurements, we deduce an average 
IRAC8.0 /M,SX8.3 flux ratio o f 0.90 ± 0.36, consistent with the median value of 1.2 ± 0.2 
fo und by ICohen et al.l (l2007bl) for PNs but slightly lower than that of 1.55 ± 0.15 obtained 
by ICohen et al.l (j2007al ) for H II regions. 



Cohen et al.l (l2007bl ) found a marginal relation between the IRAC8.0/M5 , A 7 8.3 flux ratio 



and nebular diameters which suggests that the aperture correction factor increases from 1.0 
at small nebular size (i.e. point-source) to > 1.4 at 77". This relation is not shown in 
Figure [51] because our sample does not include very extended PNs. The PNs investigated 
here have a diameter smaller than 60", and mostly smaller than 20", much more compact than 
the arcminutes-scale H II regions. Therefore, we infer that the absolute flux calibration for 
the compact sources studied here should be accurate, and we do not need to make aperture 
correction to obtain the integrated fluxes of the GLIMPSE II PNs. 

The [S IV] 10.5 /im line, which can contribute to the MSX 8.3 /im band but not to the 
IRAC 8.0 /im band, is partly responsible to the IRAC8.0/MS'X8.3 flux ratio. PNs often have 
higher excitation class than H II regions, and thus are supposed to have stronger [S IV] 
emission. As a result, the IRAC8.0/M5X8.3 flux ratio in PNs (especially high-excitation 
PNs) might be lower compared with that in H II regions. 



We cannot completely rule out variability as the partial cause of the IRAC and MSX fiux 
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difference for individual PNs since the GLIMPSE survey was performed about ten years after 
the MSX observations. If this is the case, our study will have implications to understand 
the pulsations mechanism or binary system in PNs. Another possible reason causing the 
IRAC / MSX flux discrepancy in individual PNs is that the MSX fluxes of some PNs might 
be contaminated by blended sources since MSX has a resolution approximate ten times lower 
than IRAC. 



5.3. MIPS 24 /im versus MSX 21 /im Integrated Fluxes 

The ratio of MIPS 24 /iin to MSX 21 /im fluxes provides an opportunity to establish the 
reliability of absolute diffuse calibration and our measurements of MIPS 24 /im fluxes. For 
this purpose, one should at first estimate the predicted MI PS24/M5X21 flux ratio through 



integrating the spectral response functions to the SEDs. ICohenl (120091 ) investigated the 
diffuse absolute calibration of the MIPS 24 /im channel using a sample of H II regions, 
and obtained the observed MIPS24/M5A / 21 flux ratio of 1.19 ± 0.03. Comparing with the 
predicted value of 1.18 ± 0.02, which is obtained from three SEDs pertinent to H II regions, 
they concluded that the absolute diffuse calibration of MIPS24 is consistent with that of 
MSXll within 3% (la). We reasonably suppose that predicted MIPS24/MSX21 flux ratio 
of PNs is nearly the same as that of H II regions. 

In Figure [52] we compare the MIPS 24 /im and MSX 21 /im integrated fluxes of the 
GLIMPSE II PNs. The sensitivity of the MSXll band is approximately 1000 times lower 
than of the MIPS24 band, and thus only relatively bright PNs can be detected by MSXll. 
We found that nine of the GLIMPSE II PNs have MSXll counterparts. Since the MIPS24 
band has a much lower saturated level, very bright sources in the MSX catalogue could 
be saturated in MIPS. This is clearly reflected in Figure [52j Excluding the two saturated 
sources, we obtain the MIPS24 / MS X1 1 flux ratio of 1.19 ± 0.12, in perfect agreeme nt with 



value for H II regions (ICohenll2009l ). We note that the H II regions investigated by I Cohen 
( 120091 ) are much more extended, and generally have 24 /im integrated fluxes ranging from 
10 4 -10 5 5 mJy, significantly larger than those of the GLIMPSE II PNs in Figure [521 (about 
10 3 - 2 -10 3 - 9 mJy). The consistent results obtained from PNs and H II regions, which have 
remarkably different spatial sizes and integrated fluxes, strongly validate the reliabilities of 
the absolute diffuse calibration of MIPS24 and our flux measurements. 
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5.4. Comparison with Radio Flux Densities 



A tight correlation between far- and mid- infrared and radio emission from star-forming 
galaxies h as been established a nd has been a subject of intense discussion for many years 
(see, e.g., lAppleton et al.ll2004 , and the references therein). The origin of the tightness of 



this correlation is, however, not completely understood. It would be interesting to test the 
MIR/radio relation in PNs, although they have quite different physical conditions from star 
formation regions. 

If the MIR emission is dominated by dust heated by ultraviolet (UV) photons from the 
central stars which also ionize the gaseous nebulae and are responsible for free-free thermal 
emission in the r adio continuum, we would see a correlation between the MIR and radio 
emission in PNs. I Cohen et al.l (l2007bl ) compared the IRAC 8.0 /im, MSX 8.3 /im, and radio 
flux densities and obtained an overall MIR/radio ratio of 4.6 ± 1.2. They also suggest that, 
unlike that in H II regions, the MIR/radio ratio in PNs might change with PN evolution 
as the result of decreasing amount of UV photons available to pump the AIB into MIR 
fluorescence. Even without considering the line emission factor, the dust continuum flux of 
PNs will decrease as the dust cools due to geometric dilution of the heating photons, which 
will occur at a different rate as the radio continuum flux which decreases as the electron 
density drops as the result of expansion even the electron temperature remains constant. 

We note that a ll the GLIMPSE II PNs are located within the NVSS 1.4 GHz survey 
(jCondon et al.lll998l ) field, and 23 of them have NVSS detections. They have sizes smaller 
than the 1.4 GHz beam size of VLA and are spatially unresolved at this frequency. In Fig- 
ure EH we compare the IRAC 8.0 /on and MIPS 24 /im integrated fluxes with those from 
NVSS catalogue. It is evident that albeit with a large dispersion, there is a positive corre- 
lation between the MIR and radio emission. The correlation coefficients in the logarithmic 
space are 0.66 for IRAC 8.0 /im vs NVSS 1.4 GHz, and 0.64 for MIPS 24 /im vs NVSS 
1.4 GHz. Figure [531 also shows that the number distributions of the IRAC8.0/im/1.4GHz 
and MIPS24/iin/1.4GHz flux ratios are peaked at ~ 3 and ~ 100, respectively. If we as- 
sume that these GLIMPSE II PNs are approximately at the same evolutionary stage and 
that the non-thermal emission is insignificant, the dispersion in Figure [531 should be due to 
self-absorption of radio emission. For most compact PNs, the turnover (from optically thick 
to optically thin) frequency is usually at 1-5 GHz, so it is quite possible that some of these 
nebulae are optically thick at 1.4 GHz. The object showing the largest MIR/radio flux ratio 
is PNG 359.2 + 01.2 which has a high density torus. 

A basic observable parameter qm[= log(Sm/ Siaghz), S being the flux densit y] is often 
used to characterize the MIR/radio relation. For i nstance. lAppleton et al.l (120041 ) obtained 
q 2 4 = 1.00±0.27 for galaxies with different redshifts; llbar et al.l (120081 ) found g 2 4 = 0.71±0.47 
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for a sample at high-redshifts; IWu et al.l (120081 ) obtained the ^24 values for dwarf star-forming 
galaxies generally ranging from 0.9 to 1.3 with one extreme case of 2.2. For the GLIMPSE II 
PNs, we obtain q 8 = 0.79 ± 0.52 and g 2 4 = 2.13 ± 0.37. The q S n value is consistent with 



that found for PNs in the Magellanic Clouds by iFilipovic et al. I (120091 ). The larger g 2 4 



value in PNs compared to those in star-forming galaxies is a reflection of the fact that PNs 
and galaxies have different radio and infrared emitting mechanisms. Galaxies radiate non- 
thermal radiation whereas PNs radiate thermal f-f emission. The dust components in PNs 
are heated by a single central star whereas the dust in galaxies is heated primarily by massive 
stars. 



6. CONCLUSIONS 

An investigation of the MASH PNs within the GLIMPSE II field has revealed that 
37 optically identified PNs have MIR counterparts. These PNs are distinguishable from the 
surrounding stars by their redder color and more extended angular sizes. From their infrared 
morphologies, we suggest that some of these objects (PNG 352.1-00.0, PNG 353.9+00.0, 
PNG 357.4-01.3, PNG 358.8-00.0) are unlikely to be PNs. We present the spatially integrated 
flux and size measurements of these PNs at 3.6, 4.5, 5.8, 8.0, and 24 //m. MIR color-color 
and color-magnitude diagrams are produced. Combining our results with those in other data 
archives, we construct the 0.82-100 /im SEDs for these objects. The primary findings from 
analysis of these results are summarized below. 

1. Compared to the GLIMPSE I PNs, the GLIMPSE II PNs are generally more compact 
and have lower spatially integrated fluxes, implying that these PNs towards the Galac- 
tic bulge might be generally more distant. Nevertheless, the MIR colors of GLIMPSE I 
and II PNs are very similar to each other, and are clearly different with those of stellar 
objects. 

2. Despite some GLIMPSE II PNs are invisible in the 3.6-5.8 /im bands and are only 
marginally above the IRAC 8.0 fim detection limitation, all of them are characterized 
with very strong 24 /im emission, which is likely to originate from extended cold dust 
envelopes. 

3. The nebular morphologies and fluxes detected in different wavelengths allow us to 
separate the different components in the nebulae. The 5.8 and 8.0 /im emission may 
be either more or less compact than their optical counterparts, respectively implying 
the presence of dust toruses centered on bipolar PNs or extended neutral envelopes 
around ionized gaseous nebulae. 
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4. From the anomalous colors of some objects, we infer the contributions of AIB emissions 
to the IRAC fluxes. 

5. The absolute diffuse calibrations are examined by comparing the IRAC 8.0 /im and 
MIPS 24 /iin to the MSX 8.3 and 21 /im fluxes. We obtain the IRAC8.0/MSX8.3 and 
MIPS24/MSX21 flux ratios of 0.90 ± 0.36 and 1.1 9 ±0.12, respectively. These results 



are rea sonabl y consis t ent w ith previous studies by I Cohen et al.l (j2007al ) , ICohen et al. 



fl2007bf ). and ICohenl fl2009h using data for H II regions and GLIMPSE I PNs, and 
suggest good calibrations for these slightly resolved GLIMPSE II PNs. 

We explore the MIR/radio relation for PNs. A loose correlation between the MIR 
(8.0 /im and 24 /im) and the NVSS 1.4 GHz fluxes is found, and can be statistically 
described by giR — the ratio between the flux densities at MIR and 1.4 GHz. For the 
GLIMPSE II PNs, we find g 8 . = 0.79±0.52 and q 2A = 2.13±0.37. The larger q m values 
and the larger dispersion in the MIR/radio correlation of PNs compared to those of 
star-forming galaxies are attributed to their respective different underlying radiation 
mechanisms. 



Since the bulk of GLIMPSE II PNs are located at higher Galactic latitudes (1° < \b\ < 
2°), we anticipate that a large number of PNs might have been detected in the GLIMPSE 3D 
legacy survey, which extends the latitude coverage to ±3°. The study of GLIMPSE 3D PNs 
will be presented in a subsequent paper. 
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Fig. 1.— The survey coverage of GLIMPSE I (dotted lines), GLIMPSE II (dotted-dashed 
lines), and GLIMPSE 3D (solid lines). The circles and pluses represent the MASH PNs 
inside and outside the GLIMPSE I/II survey coverage, respectively. The filled circles denote 
the MASH PNs having IRAC counterparts that are presented in this paper and Paper I. 
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Fig. 2. — Composite-color images of PNG 000.0 — 01.3 made from 3 IRAC bands. The 
3.6 /iin is shown as blue, the 5.8 /im is shown as green, and the 8.0 //m is shown as red. The 
lower panel is the same with the upper panel but overlaid with the contours of Ha emission. 
The target is put in the center of each panel. The coordinates are in units of R.A. (J2000.0 
on the horizontal scale) and decl. (J2000.0 on the vertical scale). 

Fig. 3.— As in Figure. El but for PNG 000.1 - 01.7. 

Fig. 4.— As in Figure. El but for PNG 000.3 - 01.6. 

Fig. 5.— As in Figure. El but for PNG 000.5 + 01.9. 

Fig. 6.— As in Figure. El but for PNG 000.9 - 01.0. 

Fig. 7.— As in Figure. El but for PNG 000.9 + 01.8. 

Fig. 8.— As in Figure. El but for PNG 001.0 - 01.9. 

Fig. 9.— As in Figure. El but for PNG 001.0 + 01.9. 

Fig. 10.— As in Figure. El but for PNG 001.6 - 01.1. 

Fig. 11.— As in Figure. El but for PNG 002.0 + 00.7. 

Fig. 12.— As in Figure. El but for PNG 002.1 - 00.9. 

Fig. 13.— As in Figure. El but for PNG 002.1 - 01.1. 

Fig. 14.— As in Figure. El but for PNG 002.2 + 00.5. 

Fig. 15.— As in Figure. El but for PNG 002.2 - 01.2. 

Fig. 16.— As in Figure. El but for PNG 003.4 + 01.4. 

Fig. 17.— As in Figure. El but for PNG 003.5 - 01.2. 

Fig. 18.— As in Figure. El but for PNG 003.5 + 01.3. 
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Fig. 19.— As in Figure. El but for PNG 003.6 - 01.3. 



Fig. 20.— As 



Fig. 21.— As 



Fie. 22.— As 



Fig. 23.— As 



Fig. 24.— As 



Fig. 25.— As 



Fig. 26.— As 



Fig. 27.— As 



Fig. 28.— As 



Fig. 29.— As 



Fig. 30.— As 



Fig. 31.— As 



Fig. 32.— As 



Fig. 33.— As 



Fig. 34.— As 



Fig. 35.— As 



Fig. 36.— As 



n Figure. El but for PNG 003.7 + 00.5. 



n Figure. El but for PNG 004.3 - 01.4. 



n Figure. El but for PNG 004.8 - 01.1. 



n Figure. El but for PNG 006.1 + 00.8. 



n Figure. El but for PNG 352.1 - 00.0. 



n Figure. El but for PNG 352.8 - 00.5. 



n Figure. El but for PNG 353.9 + 00.0. 



n Figure. El but for PNG 355.6 - 01.4. 



n Figure. El but for PNG 355.6 + 01.4. 



n Figure. El but for PNG 356.0 - 01.4. 



n Figure. El but for PNG 356.9 + 00.9. 



n Figure. El but for PNG 357.4 - 01.3. 



n Figure. El but for PNG 357.5 + 01.3. 



n Figure. El but for PNG 357.7 + 01.4. 



n Figure. El but for PNG 358.2 - 01.1 



n Figure. El but for PNG 358.8 - 00.0. 



n Figure. El but for PNG 359.1 - 01.7. 
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Fig. 37.— As in Figure. but for PNG 359.2 + 01.2. 



Fig. 38.— As in Figure. U but for PNG 359.3 - 00.9. 



Fig. 39.— MIPS 24/mi images of the GLIMPSE II PNs. From left to right: PNG 000.0-01.3 
and PNG 000.1 - 01.7 (1st row); PNG 000.3 - 01.6 and PNG 000.5 + 01.9 (2nd row); PNG 
000.9 - 01.0 and PNG 000.9 + 01.8 (3rd row); PNG 001.0 - 01.9 and PNG 001.0 + 01.9 (4th 
row). The target is put in the center of each panel. The coordinates are in units of R.A. 
(J2000.0 on the horizontal scale) and decl. (J2000.0 on the vertical scale). 



Fig. 40.— MIPS 24/im images of the GLIMPSE II PNs. From left to right: PNG 001.6-01.1 
and PNG 002.0 + 00.7 (1st row); PNG 002.1 - 00.9 and PNG 002.1 - 01.1 (2nd row); PNG 
002.2 + 00.5 and PNG 002.2 - 01.2 (3rd row); PNG 003.4 + 01.4 and PNG 003.5 - 01.2 (4th 
row). Details are indicated in Figure. [39l 



Fig. 41.— MIPS 24/zm images of the GLIMPSE II PNs. From left to right: PNG 003.5+01.3 
and PNG 003.6 - 01.3 (1st row); PNG 003.7 + 00.5 and PNG 004.3 - 01.4 (2nd row); PNG 
004.8 - 01.1 and PNG 006.1 + 00.8 (3rd row), PNG 352.1 - 00.0 and PNG 352.8 - 00.5 (4th 
row). Details are indicated in Figure. [391 



Fig. 42.— MIPS 24/xm images of the GLIMPSE II PNs. From left to right: PNG 353.9+00.0 
and PNG 355.6 - 01.4 (1st row); PNG 355.6 + 01.4 and PNG 356.0 - 01.4 (2nd row); PNG 
356.9 + 00.9 and PNG 357.4 - 01.3 (3rd row); PNG 357.5 + 01.3 and PNG 357.7 + 01.4 (4th 
row). Details are indicated in Figure. [391 



Fig. 43.— MIPS 24/im images of the GLIMPSE II PNs. From left to right: PNG 358.2- 
01.1 and PNG 358.8-00.0 (1st row); PNG 359.1-01.7 and PNG 359.2+01.2 (2nd row); PNG 
359.3-00.9 (3rd row). Details are indicated in Figure. [391 
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Fig. 44.— The SED of the PNs in the GLIMPSE II sample. The filled triangles, filled 
diamonds, filled circles, open circles, and filled squares are from the DENIS, 2MASS, 
GLIMPSE/MIPSGAL, IRAS, and MSX survey, respectively. In the case of PNG 359.3-00.9 
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Fig. 45. — The ISO spectrum of Hb 5. The AIB emission and strong atomic lines are 
identified. The normalized relative spectral response curves for the IRAC and MIPS 24 fxm 
bands are overlaid. 
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Fig. 46.- IRAC color-color plot ([3.6] - [4.5] vs. [ 5.8] 
GLIMPSE I and GLIMPSE II PNs. The observation results by Eora et all fl2004h and 



[8.01 ) for the 



Phillips Sz Ramos-Larios J2008al lbh and the theoretical calculations by I Cohen et al.l (]2007al ) 



are overlaid. Th e enclos ed zone by dotted lines represents that occupied by the PNs studied 
by ICohen et al.l (l2007bl ). We also show the colors of UC H II regions ( Ide la Fuente et al. 
20081 ) and three types of YSOs (star with infalling envelop e — pluses, star with a disk- 



crosses, and post T-tauri star — points; IChavarria et al.ll2008l ). The filled squares connected 
by a solid line denote the prediction of blackbody radiation with various temperatures Tgg. 
The arrow in the upper left corner denotes a reddening vector of A v = 10. 
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Fig. 47.— The [3.6] - [8.0] versus [8.0] - [24] color-color diagram for the GLIMPSE II PNs. 
Symbols are otherwise same as in Figure. |46j 
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Fig. 48. — The [3.6] versus [3.6] — [8.0] (upper panel) and [8.0] versus [3.6] — [8.0] (lower 
panel) color-magnitude diagrams for the GLIMPSE I (open circles) and GLIMPSE II (filled 
circles) PNs. 
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Fig. 49. — The histogram of the relative sizes of nebular major axises measured in different 
wavelengths. 
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Fig. 50. — The flux distribution along the major axis of the bipolar nebula PNG 359.2 + 01.2 
at Ha (dotted line) and 8.0 /xm (solid line) bands. 
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Fig. 51. — IRAC 8.0 /xm versus MSX 8.3 /xm integrated fluxes for the GLIMPSE I (open 
symbols; from Paper I) and GLIMPSE II (filled symbols) PNs. The triangles, squares, and 
circles denote the PNs with diameters (R) in the ranges of R < 10", 10" < R < 20", and 
20" < R < 60", respectively. The solid diagonal line is a y = x plot. The dotted line 
represents a linear least-squares fit. The error bar is given in the lower right corner. 
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Fig. 52.— IRAC 24 /im versus MSX 21 fim integrated fluxes for the GLIMPSE II PNs. The 
solid diagonal line is a y = x plot. The dotted line represents a linear least-squares fit. The 
error bar is given in the lower right corner. 
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Fig. 53. — Left panels: IRAC 8.0 /xm (upper panel) and MIPS 24 /im (lower panel) versus 
NVSS 1.4 GHz integrated fluxes for the GLIMPSE II PNs. The dotted lines represent a 
linear fitting. Right panels: the number distributions of the 8.0//m/1.4GHz (upper panel; 
note that the two objects with extremely large flux ratios are out of the range of the abscissa) 
and 24/zm/1.4GHz (lower panel) flux ratios. 



Table 1. Known PNs. 



Object Other Name 



PNG 000.5 + 01.9 PN K 6-7 

PNG 001.0-01.9 PNK6-35 

PNG 001.0 + 01.9 Hen 2-264, PN Kl-4, PN Sa 3-69 

PNG 002.2 + 00.5 Terz N 2337 

PNG 358.2 - 01.1 PN Al 2-L, PN Bl 1-D 

PNG 358.8 - 00.0 Terz N 2022 

PNG 359.1 - 01.7 Hen 1-191, PN M 1-29, PN Sa 2-250 

PNG 359.2 + 01.2 PN PM 1-166, 19W32 

PNG 359.3 - 00.9 Hen 2-286, PN Bl 1-E, PN Hb 5, PN Sa 2-244 



Table 2. Spitzer observations of PNs in GLIMPSE II. 





Obs. Coord. (J2000) 




Flux (mjy) 






Size (arcsec) 










Object 


R.A Decl. 


3.6/im 


4.5/^m 5.8/im 8.0/im 


24/im 


3.6/im 


4.5/j.m 5.8/tm 


8.0Atm 


24/im 


Ha" 


CQRAC 



PNG 000.0-01.3 


17 51 00.4 


-29 33 51 






9.1 


20.2 


320 






9x11 


10x13 


20 


10.0x12.0 




PNG 000.1-01.7 


17 52 49.0 


-29 41 56 








9.9 


699 








12 


20 


11.0x17.0 




PNG 000.3-01.6 


17 52 52.1 


-29 30 04 








10.4 


650 








8 


13 


5.6 




PNG 000.5+01.9 


17 39 31.4 


-27 27 47 








26.1 


660 








9x11 


15 


6.0x9.5 




PNG 000.9-01.0 


17 51 43.3 


-28 38 59 








23.1 


650 








9 


15 


6.0 




PNG 000.9+01.8 


17 40 51.7 


-27 08 48 








12.0 


153 








4x6 


11 


5.5x7.0 




PNG 001.0-01.9 


17 55 43.9 


-29 04 05 








16.9 


1380 








11x15 


35 


15.0x17.0 




PNG 001.0+01.9 


17 40 26.8 


-27 01 09 










230 










13 


53.5x78.2 




PNG 001.6-01.1 


17 54 13.4 


-28 05 17 








30.3 


524 








9x10 


27 


11.0x15.0 




PNG 002.0+00.7 


17 47 28.3 


-26 49 48 


10.1 


14.4 


21.6 


39.1 


970 


4 


4 


7 


8 


25 


4.0 


0.69 


PNG 002.1-00.9 


17 54 18.8 


-27 36 36 


48.7 


62.6 


87.9 


216.2 


930 


6 


6 


9 


12 


28 


6.0 


0.85 


PNG 002.1-01.1 


17 55 10.4 


-27 41 40 






7.3 


14.7 


660 






4 


6 


31 


8.0 




PNG 002.2+00.5 


17 48 45.6 


-26 43 31 


21.0 


42.0 


102.6 


239.0 


2295 






23x32 


26x35 


32 


17.0x28.0 


2.08 


PNG 002.2-01.2 


17 55 45.6 


-27 39 41 




5.4 


5.8 


22.3 


530 




4x5 


5x7 


10x12 


32 


11.0x14.5 




PNG 003.4+01.4 


17 48 15.6 


-25 15 14 








19.4 


67 








10 


14 


12.5x14.0 




PNG 003.5-01.2 


17 58 37.0 


-26 28 47 


78.0 


108.8 


339.8 


1080.1 


2510 


6 


6 


11 


11 


30 


4.0x5.0 


2.44 


PNG 003.5+01.3 


17 48 41.7 


-25 11 34 


15.0 


17.1 


16.2 


52.4 


2260 


5 


5 


9 


10 


30 


4.0x5.0 


0.47 


PNG 003.6-01.3 


17 59 11.0 


-26 30 23 


13.2 


17.5 


15.0 


23.1 


1310 


6 


7 


8 


10 


30 


14.0 


-0.41 


PNG 003.7+00. 5^ 


17 52 08.9 


-25 27 43 










< 1660 












10.2 




PNG 004.3-01.4 c 


18 01 18.9 


-25 53 21 


3.1 


7.0 


4.4 


29.2 


1680 


4 


4 


5 


8 


34 


5.0x6.0 


1.44 


PNG 004.8-01.1 


18 01 17.0 


-25 22 37 


20.6 


31.4 


25.3 


120.6 


6205 


5x8 


5x8 


7x10 


10x13 


35 


6.0x8.0 


1.00 


PNG 006.1+00.8 


17 56 33.2 


-23 11 47 


13.8 


11.4 


7.2 


26.6 


510 








11 


32 


9.0 


-0.32 


PNG 352.1-00.0 d > e 


17 25 33.4 


-35 36 25 


678 


1412 


4873 


10599 


23000 






73x96 


107x156 


40x50 


45.0 


2.55 


PNG 352.8-00.5 / 


17 29 37.5 


-35 13 46 


57.9 


58.8 


57.5 


441.9 


> 10580 


5x7 


7x10 


13x17 


13x17 


34 


9.0x18.0 


1.43 


PNG 353.9+00.0° 


17 30 29.5 


-34 00 37 


54.0 


52.9 


50.1 


58.7 


90 


5 


6 


7 


10 


11 


4.0x5.0 


-0.91 


PNG 355.6-01.4 


17 40 55.2 


-33 24 18 








11.6 


360 








6x9 


20 


9.4 




PNG 355.6+01.4 


17 29 11.0 


-31 52 45 


9.8 


14.8 


23.5 


85.5 


2300 


4 


4 


7 


10 


30 


5.0x6.0 


1.67 


PNG 356.0-01.4 


17 41 33.4 


-33 02 15 






7.2 


28.9 


680 






5 


11 


28 


7.0x8.0 




PNG 356.9+00.9 


17 34 34.3 


-31 02 08 


16.5 


19.8 


55.8 


206.8 


4400 


5 


5 


10 


10 


30 


5.0 


2.31 


PNG 357.4-01.3 e - 9 


17 44 56.2 


-31 49 19 
















300 


300 


300 


256 




PNG 357.5+01.3 


17 34 26.0 


-30 15 18 


7.7 


6.4 


33.0 


94.7 


390 


4 


4 


8 


10 


15 


4.0x7.0 


2.51 


PNG 357.7+01.4 


17 34 46.6 


-30 04 21 


1.8 


2.2 


3.5 


10.9 


180 






4 


6 


15 


6.0x10.0 


1.27 


PNG 358.2-01.1 


17 46 03.1 


-31 03 36 


28.2 


30.6 


36.7 


55.7 


400 






17 


19 


25 


18.8x22.6 


-0.14 



Table 2 — Continued 



Object 


Obs. Coord. (J2000) 






Flux (mjy) 








Size (; 


arcsec) 








«IRAC 


R.A 


Decl. 


3.6/mi 


4.5/^m 


5.8/xm 


8.0^m 


24/xm 


3. 6^tm 


4.5/tm 


5.8/xm 


8.0£tm 


24/im 




PNG 358.8-00.0 e ' f 


17 42 42.6 


-29 51 35 


2823 


2582 


18057 


50662 


> 37000 


50 


50 


70 


80 


55 


11.3x28.3 


3.04 


PNG 359.1-01.7 


17 50 18.0 


-30 34 55 


36.4 


69.5 


63.3 


275.3 


6130 


9 


9 


13 


15 


35 


17.0x18.8 


1.31 


PNG 359.2+01.2 


17 39 02.9 


-28 56 37 


442.1 


732.0 


1361.6 


2119.6 


8450 


8x16 


8x16 


9x17 


13x22 


30x42 


9.5x26.5 


0.99 


PNG 359.3-00.9-f>' 1 


17 47 56.1 


-29 59 41 


303.1 


312.6 


354.8 


> 6369 


> 18800 


12 


12 


23 


23 


30 


28.3x60.3 





a Taken from lParker et alj l l200ll. l2006ri and lMiszalaki et al. I <]2008h . 

b Badly blended with a bright star, which saturates all the IRAC bands. 

c The MIPS image is probably blended with two unresolved field stars. 

d The Spitzer images reveal diffuse tail. The flux and size of the central bright region are given for the MIPS24 band. 
c Non-PNs. 

'The MIPS images are saturated. 

s It shows large circle structure with ambiguous edge, and the given sizes are rough. 
h The 8.0 £im image is saturated. 



Table 3. Other flux measurements 



DENIS 2MASS MSX IRAS 





I 


J 


K 


J 


H 


K 


8.28 (im 


12.13 urn 


14.65 Atm 


21.3 urn 


12 jim 


25 fj,m 


60 fim 


100 Atm 


Object 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


PNG 002.0+00.7 




13.902 


11.764 


13.858 


12.486 


11.82 


















PNG 002.1-00.9 


15.616 


14.268 


12.009 


13.827 


12.702 


11.56 


0.159 




1.121 












PNG 003.5—01.2 


16.089 


12.894 


10.980 


13.020 


12.264 


11.03 


0.998 


2.653 


2.153 


2.266 


1.996 


3.438 


9.337 


43.230 


PNG 003.5+01.3 














0.195 




1.450 


2.052 










PNG 003.6-01.3 


15.963 


13.137 


11.606 


13.123 


11.978 


11.59 


















PNG 004.3-01.4 






















2.272 


1.864 


2.691 


3.623 


PNG 004.8-01.1 


17.951 


14.203 


12.166 


14.333 


12.815 


11.71 


0.258 




3.162 


5.381 


3.329 


8.839 


20.550 


106.900 


PNG 352.8-00.5 


18.071 


14.654 


12.250 


14.995 


13.929 


12.24 


0.838 


1.774 


4.84 


13.878 


1.693 


25.430 


137.400 


172.900 


PNG 353.9+00.0 








13.999 


12.284 


10.97 


















PNG 355.6+01.4 




14.002 


12.361 


13.865 


12.708 


12.04 


0.113 




0.725 


1.720 










PNG 356.9+00.9 








14.495 


13.384 


12.34 


0.181 


1.946 


0.635 


4.195 










PNG 358.8-00.0 








10.375 


9.763 


9.285 


















PNG 359.1-01.7 














0.532 




5.100 


4.500 










PNG 359.2+01.2 








11.442 


10.283 


9.550 


1.908 


3.452 


4.459 


7.076 










PNG 359.3-00.9 


13.149 


11.265 


9.980 


11.037 


10.852 


9.854 


4.985 


11.118 


28.664 


48.165 


11.680 


79.240 


134.500 


311.800 



